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Abstract 

We consider the new-physics (NP) solution to the polarization puzzle m B ^ (pK* 
decays. We note that any such solution must reproduce the data in 4>Ks, where 

there are disagreements with the standard model in CP-asymmetry measurements. 
We examine 10 NP operators, of S/P, V/A and T variety. We find that, as long 
as 5° —>■ (j)Ks exhibits large CP-violating effects, no single operator can explain the 
observations in both B (j)K* and B^ <pKs. For 2-NP-operator solutions, there 
are four possibilities, all of S/P type, which are presently allowed. We discuss ways 
of distinguishing among these solutions in the future. Models which contain only V/A 
operators, such as those with supersymmetry or extra Z' bosons, cannot explain both 
B^ — >■ and B — > (pK* data. On the other hand, the two-Higgs-doublet model, 
which has only S/P operators, is favored. 
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5 — * V1V2 decays (the Vi are vector mesons) really represent three transitions when 
the spins of the Vi, are taken into account. That is, the final state can be transversely 
(2 states) or longitudinally (1 state) polarized. A naive calculation within the standard 
model (SM) shows that the transverse amplitudes are suppressed by a factor of size 
mv/rriB {V is one of the vector mesons) with respect to the longitudinal amplitude. 
One then expects the fraction of transverse decays, Jt, to be much less than the fraction 
of longitudinal decays, f^- 

However, it was observed that these two fractions are roughly equal in the decay 
B (j)K*: fr/fL — 1 [HEIE]- This is the "polarization puzzle." If one goes beyond 
the naive SM, there are two explanations [3] which account for this surprising result: 
penguin annihilation (PA) [5] and non-perturbative rescattering [6l[7]. Still, there are 
question marks associated with both of these. First, PA is a subleading amplitude that 
is power suppressed by 0(1/ nib). Second, for rescattering, it is not obvious whether 
such a non-perturbative effect is of leading or subleading order. Hence, the SM expla- 
nations of the large /t//l generally require enhanced subleading amplitudes (certain 
for PA; possible for rescattering). 

One can also explain the /t//l measurement by introducing physics beyond the 
SM. Suppose there are new-physics (NP) contributions to the b sss quark-level 
amplitude. If their form is chosen correctly, these will contribute dominantly to in 
B — > (pK* and not to f^, so that one can reproduce the measured value of /t//l if the 
NP amplitude has the right size [8]. In this paper, we explore the NP explanation 
- we assume that neither PA nor non-perturbative rescattering produce dominant 
contributions to the transverse amplitudes, and are therefore not the explanation of 
the measurement of a large /t//l- 



BR 


8.3+}-^ X 10-^ 






-A-CP 


0.23 ±0.15 



Table 1: Measurements of B^ <pKs- Included are the branching fraction (BR) 
[Ulinilll], the indirect (mixing-induced) CP asymmetry {Sep) [HI [12], and the direct 
CP symmetry (A^p) W^T^- 

Now, any NP contribution to 6 — sss will also affect B^ <pKs. Thus, any 
constraints on such NP must take into account the measurements of both B (pK* and 
B^ — i> (pKg. The B'^ <pKs data are shown in Table 1. The CP- violating observables 
are particularly intriguing. Within the SM, including small corrections, the indirect 
(mixing-induced) CP asymmetry in B^ (pKg {Sqp) is expected to be a bit larger 
than that in charmonium B^ decays [T3], found to be 5'cp(charmonium) = 0.672±0.024 
pT] . In addition, the direct CP asymmetry in B^ — »• (pKg {Acp) is expected to vanish. 
In other words, the central values of both of these measurements exhibit disagreements 
with the expectations of the SM. This provides a hint of NP in 6 — > s transition^. 

"^A more significant hint (signal?) of NP, also in 6 — > s transitions, is provided by i? — > ttK decays, 
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On the other hand, the errors are sufficiently large that the discrepancies are only at 
the level of < 2a . This means that the hint of NP is not statistically significant. It also 
means that any constraints on NP in i? — > 07^^* are not that strong. For this reason, 
in this paper we also perform the analysis in the scenario in which the hint of NP in 

— > (^Ks becomes a true signal in the future. That is, in this case we assume that 
the future measurements of Sep and Acp stay at their present central values, but the 
errors are reduced by a factor of 2. If the discrepancies with the SM in i?^ — »• ^K^ get 
more pronounced, it will be necessary to consider the results from this second scenario. 

We assume that there is a NP contribution to 6 ^ sss. We consider operators of 
the form 

orf = ^7^(1 - 75)6 S7'^(i-75)s, 

Or^" = S7^(l-75)&^7''(l + 75)s, 

Oil" = S7,(l + 75)&^7'^(l-75)s, 

0];^f = S7,(l + 75)&^7'^(l + 75)s, 

Off = s(l-75)&^-(l-75>, 

Of/ = ^(1-75)6^(1 + 75)5, 

Off = ^(1 + 75)6^(1-75)5, 

Off = s(l + 75)&5fl+75)s, 
01 = sa^,(l-75)6sa^^(l-75)s, 

01 = sa^,(l + 75)6 sa^^(l + 75)5. (1) 

In the above operators, we take the colors of the quark fields in each current to be the 
same. This is the case in most typical NP models (multi-Higgs-doublets, supersymme- 
try, extra Z's, etc.). For S/P operators, a Fierz transformation of the fermions and 
colors is required in order to get a non- vanishing contribution to the production of the 
final-state vector meson from the vacuum (within factorization): 

Off,,) = -^s-(l±75)fesfl±75)s-^s-a,.(l±75)&5-^^'^(l±75>, 
Off,«, = -^^7m(1±75)& 5-7^(1 + 75)5. (2) 

(We have neglected the octet piece coming from the color Fierz transformation, which 
is justified within factorization.) 

There are also 10 operators in which one has different quark colors in the currents. 
However, these can be obtained from Eq. ([T]) as follows. Suppose that there is only one 
type of Lorentz structure. The effective Hamiltonian can then be written 

Heff = BisOibs02S + B2SaOib^S(302Sa , (3) 



where the disagreement with the SM has reached the 5a level jTTj, assuming that \C/T\ is small, as 
is expected in the SM [Hj (C and T are diagrams contributing to the decay [T5]). 



2 



where Bi^2 are complex coefficients, Oi^2 represent the Lorentz structure [S/P, V/A or 
T), and a,f5 are color indices. We will call the operator associated with Bi "color al- 
lowed" and that associated with B2 "color suppressed." In factorization, the amplitude 
for B (f)K [K = Kg or K*) has the following structure: 

A{B ^ <i)K) = ^1 + ^2, (4) 

where 

A, = {B, + B2/N,){K\-sOrh\B){(l>\-s02s\()), 

A2 = iB2 + Bi/N,){K\sO,^b\B){^\sO,^s\0) . (5) 

In the above, sOipbsOipS is obtained from sOibs02S by performing a Fierz transfor- 
mation of the fermions and the colors. The color octet piece is neglected, because it 
does not lead to the production of a from the vacuum. 

Now, if Ai or A2 vanishes, or if s0^phs02FS is the same as sOibs02S, then there is 
only one amplitude, and we can work only with color-allowed operators with a general 
coefficient - the color-suppressed operators are implicitly included in them. As we 
show below, this holds for all the operators of Eq. ([1]). In what follows, the key point 
is that any S/P operator does not contribute to the decay because it cannot produce 
a from the vacuum. However, V/A and T operators do give a nonzero contribution. 

• Lorentz structure (V ± A) x (V ± A): sO^phsOipS = sOibs02S. 

• {V ± A) X [V ^ A): Fierz transforms into an S/P operator. Thus, A2 = 0. 

• (5 ± P) X (S* =F P): Fierz transforms into a V/A operator. Thus, Ai = 0. 

• (5* ± P) X (S* ± P): Fierz transforms into a combination of an S/P and a T 
operator. Thus, Ai = and A2 7^ 0. 

• T: Fierz transforms into a combination of an S/P and a T operator. Thus, we 
effectively have sOipbsOipS = sOibs02S. 

In all cases, there is only one amplitude in Eq. ([5]) above, and so the operators of 
Eq. ([1]) contain all the "color-suppressed" operators. 

We begin by examining the case where a single NP operator is added, contributing 
to the b — >■ sss amplitude. As noted above, this affects both B (pK* and B^ 4>Ks, 
and we compute the order of magnitude of the contribution of each of the NP operators 
to these decays as follows. Consider 0^{^ (for the orders of magnitude, we ignore the 
1/Nc coming from the inclusion of the color-suppressed operators): 

(</)X|0,T|p°) = {<f>K\-s^,{l-^,)b-sY{l-l5)s\B', 

^ {<j>K\{V-A)®{V-A)\B',) 

= {K\V-A\b',) {<j>\V-A\or 

= {K\ V P°), (01 V \0r - {K\ A Bf)^ (01 V |0)^ . (6) 
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We can now use the factorized matrix elements discussed in appendix El We obtain 



{<pKs\0 



V/A 
LL 



{(f)K*\0 
{(f)K*\0 



V/A 
LL 



V/A 
LL 



B', 



B 



A=0 



A=+ 



A=- 



2/+/^m,p, = 0(l), 



2m 



K* 



I I 2 \ A 



Am\p'lA2 



-2if^V— — h if^m^[mB + mK*>Ai 



itlb + niK* 

2\ 



o((AqcdM) 

+2if^V- 



rriB + rriK* 
0{AQCD/mb) , 



+ if4>m^{mB + mK,)Ai 



(7) 



where Pc is the magnitude of the momentum of final-state particles in the B'^ rest frame. 
The values of /+, Ai, A2 and V are given in appendix |X1 Note: the polarizations 
in B ^ (f)K* are denoted L (longitudinal) and ||, _L (transverse). However, above 
we refer to + and — polarizations - the transverse (^|j,±) and helicity amplitudes 
(A±) are related by A^^^ = {A+ ± A^) /V2 (and Ay,^ = (A_ ± A+) /V2 for the CP- 
conjugate amplitudes). The above B^ (pK amplitudes correspond roughly to orders 
of magnitude 1, 1, ^, where ^ = 0{AQCD/fnb)- 







0i^*(L) 






r)v/A r)v/A 

^ LL ) ^ LR 

r)V/A rtv/A 

^ RL 1 ^ RR 
QS/P 
^ LL 
QS/P 
^ LR 
QS/P 
^ RL 
QS/P 
^ RR 

01 

QT 
^ R 


1 
1 

e 
1 

1 
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1 
1 

1 
1 


e 
1 

e 

e 
1 

e 


e 
e 

e 
1 

e 
1 



Table 2: Relative orders of magnitude of the contribution of NP operators to the 
amplitudes of B^ — > (pKs and the three polarizations of -B — > (pK* = O^Agcn/^b)]- 

We have analyzed all NP operators similarly. The results are shown in Table [2j 
In order to generate large CP asymmetries in B^ (p^s, the contribution of the NP 
operator must be large [0(1)]. This points to the four O^^'* operators, Of^^ or O^'^^. 
In order to reproduce the /t//l measurement in B —>■ (pK*, the NP contribution to 
a transverse polarization must be large [0(1)], while not contributing significantly to 
the longitudinal polarization. We see that only NP operators of the form Of^^, O^'^, 
01 or Ob satisfy this criterion [8]. 



4 



However, the key point here is that there is no NP operator that significantly affects 
both B'^ — > (pKg and a transverse amphtude of i? — >■ (j)K*. We therefore conclude that, 
as long as —>■ (pKs exhibits large CP-violating effects, there is no single NP operator 
which can account for the observations in both B'^ — > (pKg and B (pK* decays. Note: 
if one assumes that there are no NP signals in this decay, which might be justified with 
the present errors, then 1-NP-operator solutions are still possible. But if one assumes 
that there are NP signals here, as would clearly be indicated if the errors are reduced 
by a factor of 2, then 1-NP-operator solutions are not possible. 

We now turn to the case where two NP operators are added. Here it is obvious 
that one of the operators must contribute significantly to B^ —>■ (pK^ (6 possibilities), 
and the other to a transverse amplitude of S — > (pK* (4 possibilities). Thus, we must 
in principle consider 24 pairs of operators. However, this number can be reduced as 
follows. In any reasonable NP model, if there are only two new operators, these are 
typically both of the V/A, S/P or T variety. There are no pairs of V/A or T operators 
which give large contributions to both decays, so there are only 4 pairs of NP operators 
to examine: (Off, Of//), (Of/, Off), (Of//, Of/) and (Of/, Of/). The results of 
Table [2] give only the general size of contributions, so it is necessary to perform a fit to 
see which pairs of NP operators can account for the observations in both decays, and 
to what extent. We do this below. 



BR 


9.8^^;^ X 10-^ 


Acp 


0.01 ±0.05 


h 


0.480 ± 0.030 


h 


0.241 ±0.029 



Table 3: Measurements of B ^ cpK*^ . Included are the branching fraction (BR), the 
direct CP symmetry (A^p), and the fraction of longitudinal and ± decays, and /_l 

The fit includes the three observables of B^ —>■ (pKg shown in Table [TJ BR, Sep, 
Ac p. It also includes four observables of B ^ cpK*: BR, Acp, /l, f±- The latest values 
are given in Table 1^. (There are other measurements of i? ^ 4>K*, but they are not 
used in this paper.) 

The NP contributions to the above quantities are taken into account as follows. We 
are considering the effect of the SM and two NP operators (we generally refer to them as 
Ol and O2). The strength of the NP is parametrized by unknown complex coefficients 
(referred to as Ci and O2). The SM piece is calculable within QCD factorization 
(QCDf) [17], and we take the value of its contribution from there. Ci and C2 generally 
each have a magnitude, a weak phase, and a strong phase. However, in Ref. [18], it was 
shown that the NP strong phases are negligible compared to that of the (dominant) 

^In Table El it is given that /z, = 0.48, f± = 0.24. The /, are defined such that /|| = 1 - - = 
0.28, so that indeed U/f, = (/^ + /n)//^ ~ 1. 



5 



SM contribution. Thus, we have only four free parameters: two NP magnitudes and 
two NP weak phases of Ci and C2. The decay amphtudes for a given hehcity A can 
be written in terms of these free parameters simply by computing the matrix elements 
for —>■ (pKs and each polarization of -B ^ (j)K*. In general, they take the form 



= A^,,, + {(f)K\CiOi + C202 



2 5 



where 2 are the factorized matrix elements given in appendix |Al The CP-conjugate 
amplitudes A are obtained by changing the sign of the unknown weak phases in Ci 
and C2. 

All observables can be expressed in terms of the amplitudes A^. The branching 
ratio is given by 

BR=-^^y\A^W (9) 

where Tg is the lifetime of the meson. This applies to both B^ — > (l)Kg and B 
(f)K* . The time- independent CP asymmetry is given by 

For B^ 4'Ks, A is the amplitude and A is its CP-conjugate. For B (pK*, 
|y4p = J2x lA^l"^. The time-dependent CP-asymmetry in B'^ (pKg is given by 

_ ^ lmie-'^^A*A) 
^^---^ + • ^^^^ 

Finally, the helicity fractions oi B ^ (j)K* are defined as usual by 

|y^i,ll,X|2 

and similarly for the CP-conjugates n x- Combining them, fL,\\,± are defined by 



/.,ii,x = ^(/."„,, + /f,„J. (13) 



The 7 observables can therefore be expressed in terms of the 4 free parameters 
of Ci and C2. Thus, it is possible to perform a fit to obtain the preferred values 
of these parameters, and to determine whether it is possible to account for the data 
with the addition of certain NP operators. However, there is a complication in all of 
this. In the SM, the decays B^ (pKg and B — ^ (pK* are dominated by the b —>■ s 
penguin amplitude, P'. P' is actually composed of three pieces, P^, P^ and P/, where 
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the subscript refers to the internal quark in the loop (the pieces P^^ are rescattering 
amplitudes generated mainly from tree- level operators). We can write 

p' = v:,v^,Pi + v:,v,,p', + v:,VtsP[ 
^ (P.' - p;) . (14) 

In writing the second line, we have used the unitarity of the Cabibbo-Kobayashi- 
Maskawa matrix to ehminate the V^lVts term, and we have dropped the V*i,Vus term 
since ^ iKt^sl- As the weak phase of V*fjVcs is zero, P' has only its magnitude 

and a strong phase. It is these quantities that are calculated in QCDf. Unfortunately, 
as we see below, the QCDf results are not very precise. 

The QCDf calculation is treated by applying Refs. [TTl [191 120! straightforwardly. 
For (pKs 



-Si 



P' — P' 1 1 

where the A^j^^ are form factors as defined in appendix [Aj and the a's are defined in 
Ref. [19]. (Note: the above formula could contain (5 terms. However, we have neglected 
all such pieces, consistent with our assumption that PA is not present.) For B 
the SM penguin has the same form, but with explicit polarization dependence A: 

p'A _ p'A 11 

-±- = a%'^ + a^^ - -d"'^ - -a^'^ = , fl6) 

The a'^'s are defined in Ref. |[20j. Several inputs are required in order to get magnitudes 
and strong phases. For quark masses, BBNS values [T7] were used, allowing them to 
vary within a range of la. For meson masses. Particle Data Group values |9j were used, 
fixed at their central value. Wilson coefficients were calculated using Refs. [211 [22] . with 
the renormalization scale /z allowed to vary within [m{,/2,mb]. For decay constants, 
values from Table [TUl in appendix \K\ were used, within a range of la. For form factors, 
fixed values of Table [HI in appendix [X] were used, but we studied each of the three 
cases (minimal, central and maximal values). The allowed ranges of the SM penguin 
amplitudes are summarized in Table [H The SM penguin of negative helicity is neglected 
because of the small form factors (V^j^* = V^j^* ~ 0). As expected, form factors have 
little impact on the values of the SM penguin amplitudes V^k^ 'P^k*i since they 
contribute at subleading order. Numerical variations are mainly due to the random 
scanning of the parameter space. 

Ideally, in order to take into account the ranges of the QCDf determinations, one 
would scan over the allowed regions of all magnitudes and strong phases. For each set 
of SM values, the would be evaluated. In this way, we could find the best fit (i.e. 
smallest value of xLin) ^ach of the 2-NP-operator solutions. Unfortunately, this is 
not possible, as the space of SM values is too large (e.g. if we take 10 SM values/region. 
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Minimal 


Central 


Maximal 




[0.031,0.064] 


[0.031,0.062] 


[0.031,0.060] 


Arg(^0Xs)(rad) 


[3.2,3.6] 


[3.2,3.6] 


[3.2,3.6] 


1' <j,K*\ 


[0.025,0.036] 


[0.026,0.036] 


[0.027,0.037] 


Arg(po^O(rad) 


[3.4,3.6] 


[3.4,3.6] 


[3.4,3.6] 




[0.031,0.062] 


[0.033,0.062] 


[0.033,0.061] 


Arg(P+,0(rad) 


[3.0,3.4] 


[3.0,3.4] 


[3.0,3.4] 



Table 4: Allowed ranges for SM penguin magnitudes and strong phases according to 
QCDf, for the three sets of form-factor values. 

we would require 10® evaluations). As a compromise, we have adopted the following 
procedure: we fix the SM strong phases to their central values, and scan over the SM 
magnitudes. However, we have checked what happens when we take different values 
for the strong phases. We find that the numbers can change quite a bit, but a bad 
fit cannot be turned into a good fit. 



Operators 


Minimal 


Central 


Maximal 


(Of//, Of//) 
(Of/, Of/) 
(Of/, Of/) 
(Of/, Of/) 

LL{RR)1 '^RL(LR)) 

(01,01) 


2.6 (45.7%) 
1.4 (70.6%) 
1.9 (59.3%) 

1.7 (63.7%) 
15.7 (0.13%) 
3.6 (30.8%) 


2.8 (42.4%) 
1.3 (72.9%) 
1.7 (63.7%) 
1.7 (63.7%) 
10.6 (1.4%) 
3.6 (30.8%) 


3.1 (37.6%) 
1.3 (72.9%) 
1.6 (65.9%) 
1.6 (65.9%) 
7.1 (6.9%) 
3.9 (27.2%) 



Table 5: Best-fit results (Xmm) ^^r pairs of NP operators, with present-day errors on 
Sep and Acp in (pKg. The calculation was done for the three sets of form factors 

(minimal, central and maximal). 

The results are shown in Tables [5] (current errors on Sep and Acp) and [6] (errors on 
Sep and Acp reduced by a factor of 2). Here we present the smallest value of Xmin (best 
fit) for each of the 2-NP-operator solutions. (The number in parentheses indicates the 
quality of the fit, and depends on Xmin and d.o.f. individually. 50% or more is a good 
fit; fits which are substantially less than 50% are poorer.) In all cases, the worst fit is 
given by a large value of Xmin^ with a 0% quahty of fit. 

From Table O we see that, with current errors on Sep and Acp, all four S/P 2-NP- 
operator solutions give good fits to the 5° — >• (pKg and B — > (j)K* data. We also show 
the fit results for V/A and T 2-NP-operator solutions. We see that the V/A solution 
gives a very poor fit, but the T solution, while not as good as any of the S/P pairs, is 
still acceptable. 

If the errors on Scp and Acp are reduced by a factor of 2 (Table [6]), we find that 
no S/P 2-NP-operator hypothesis is an excellent fit to the data. On the other hand. 
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Operators 


Minimal 


Central 


Maximal 


(01,01) 


6.3 (9.8%) 
4.3 (23.1%) 
5.2 (15.8%) 
4.9 (17.9%) 
20.3 (0.01%) 
13.7 (0.33%) 


7.4 (6.0%) 
4.0 (26.1%) 
5.8 (12.2%) 
4.7 (19.5%) 
15.9 (0.12%) 
13.5 (0.37%) 


8.6 (3.5%) 
3.9 (27.2%) 
5.6 (13.3%) 
4.5 (21.2%) 
10.9 (1.2%) 
14.0 (0.29%) 



Table 6: Best-fit results (Xmin) pairs of NP operators, in the scenario in which the 
errors on Scp and Acp in (pKs are reduced by a factor of 2. The calculation was 

done for the three sets of form factors (minimal, central and maximal). 

none of them is ruled out, either. The most that one can say is that (Of^^, O^^^) and 
(Orr, Oy^^) are disfavored, but even this is not very strong. On the other hand, in 
this case both V/A and T 2-NP-operator solutions are essentially ruled out. 

In both error scenarios, it is the large direct CP-asymmetry measurement in — >■ 
(pKg which is hardest to accommodate. It will be important to pay attention to this 
observable in the future to determine which NP solutions are viable. 

The fact that the best fit and worst fit have substantially different Xmm shows that 
the contribution from the SM is significant, and that all ranges would be reduced 
quite a bit with an improved determination of the SM values. In fact, one could easily 
obtain poor fits for all pairs of NP operators (as well as the SM). 

Above, we have shown that all 2-NP-operator solutions involving S/P operators 
are viable, but those which contain only V/A ot T operators are disfavored or ruled 
out. Obviously, any realistic NP model which contains more than 2 operators - and 
most do - will also be allowed, as long as the observations in both B^ —>■ (pKg and 
B —>■ (f)K* decays are explained. In Table [3, we show which types of operators are 
present for some simple NP models. Even though models with supersymmetrjf^ or 
extra Z' bosons typically generate several operators, they are all of V/A type. As such, 
they cannot explain both B^ 4>Ks and B (pK* data. On the other hand, the 
two-Higgs-doublet model is favored because it contains only S/P operators (perhaps 
all 4 pairs), and can potentially accomodate both B^ <pKs and B (pK*. 



Models 


V/A 


S/P 


T 


Supersymmetry [23l 


[21 [25] 


X 






Two Higgs doublets 


[261 [27] 




X 




Extra Z' bosons 


[2H] 


X 







Table 7: Summary of operator content for some simple NP models. 

Even with the assumption of reduced errors on the CP-violating observables in 
^Models with supersymmetry generate the b — > sss transition mainly through squark-gluino loops. 

9 



— ^ 4>Ks, all four S/P 2-NP-operator solutions are allowed. This raises the obvious 
question: is there any way to distinguish these solutions? Clearly, smaller errors on the 
experimental measurements and/or the theoretical determination of the SM contribu- 
tion can help. With these, it may be that the Xmm of solution is strongly preferred 
over that of the other three. 

However, there is an additional possibility. In any B V1V2 decay, one can con- 
struct the triple product (TP). In the rest frame of the B, the TP takes the form 
q-{6iX 62), where g is the momentum of one of the final vector mesons, and Ei and £2 
are the polarizations of Vi and V2. There are two TP's, which can be written 



lm{A^A*) 



|Ao|2 + |A|||2+|/lj 



Im{A±A*, 



|Ao|2 + |A|,|2+|AxP 



(17) 



The At''^^ for B decays are defined similarly. The TP asymmetry is defined 



Aii) 



The "fake" TP asymmetry At^ is given by the same definition, but the TP's are added 
rather than subtracted. Note that the fake TP asymmetry can be nonzero even if CP 
is conserved. 



Operators 










(ofr,o^//) 
(ofr,of/f) 


[-0.30, -0.27] 
[0.29, 0.32] 
[0.26, 0.28] 

[-0.33, -0.31] 


[0.030, 0.062] 
[-0.008, 0.014] 
[-0.099, 0.056] 
[-0.036, -0.011] 


[0.16, 0.22] 
[-0.17, -0.14] 
[-0.037, 0.090] 
[-0.001, 0.000] 


[-0.006, -0.004] 
[-0.003, 0.000] 
[-0.004, 0.001] 
0.000 



Table 8: Predictions of all S/P 2-NP-operator solutions of Table [5] for the central 
values of the real and fake TP asymmetries in B ^ (pK* . The ranges of TP-asymmetry 
predictions correspond to the full variation of form-factor values. 



The above applies to B (f)K*. In Table |8] we present the central values of At 
and At ' , calculated for each of the S/ P solutions shown in Table [51 The ranges 
correspond to the sets of form-factor values varying from minimal to maximal. We do 
not include errors because they are very large with the current data. In any case, our 
point in presenting the results of Tabled is the following. It is clear that different S/P 
2-NP-operator solutions lead to very different patterns of central values of predictions 
for the TP asymmetries. This emphasizes the usefulness of TP's for distinguishing 

^Note: in contrast to Ref. [52], this definition involves a subtraction rather than an addition. This 
is because we have defined and in such a way that At'' is zero in the absence of CP 

violation. 
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the various NP solutions, and we strongly encourage experimentalists to make such 
measurements. 

In summary, a "polarization puzzle" has been observed in S — > 0^*, namely that 
the fraction of transversely-polarized decays is about equal to that of longitudinally- 
polarized decays, in contrast to the expectations of the naive standard model (SM). 
In this paper, we explore the new-physics (NP) solution to this puzzle. We first note 
that any NP explanation must also be consistent with the observations in — > (pK^. 
This decay is particularly intriguing since the present measurements of CP- violating 
asymmetries are in disagreement with the SM. On the other hand, the errors are 
still sufficiently large that this discrepancy is not statistically significant. As such, 
any constraints on NP in 5 — > (pK* are not that stringent. For this reason, we 
also perform the analysis with the assumption that future measurements will show a 
greater statistical discrepancy. That is, we use the central values of the CP-asymmetry 
measurements, but take the errors to be reduced by a factor of 2. 

We consider 10 NP operators, of types S/P, V/A and T. We first show that, as long 
as 5^ <pKs exhibits large CP- violating effects, no single NP operator can explain the 
data in both B — > (pK* and B^ — > (pKg decays. Turning to 2-NP-operator solutions, 
it is clear that one of the operators must contribute significantly to B^ (pKg, and 
the other to a transverse amplitude of B ^ (f)K*. In any realistic NP model the two 
operators are typically both of the V/A, S/P or T type. However, no pairs of V/A or T 
operators give large contributions to both decays, so that only 4 pairs of S/ P operators 
need be considered. We have performed fits to several observables in S — > (pK* and 
4>Ks decays, and find that all four 2-NP-operator solutions are allowed. (We 
also show explicitly that the V/A and T solutions are disfavored or ruled out.) 

One can distinguish among the solutions in several ways. If the experimental errors 
on future measurements are improved, one solution might be preferred. Alternatively, 
the theoretical uncertainty can be reduced if the SM contribution to B ^ (j)K* and 
i?2 <pKg is better determined. Finally, the four solutions predict a very different 
pattern of triple-product asymmetries in i? — > (pK*. Their measurement could help 
distinguish among the possible NP solutions. 

Finally, any realistic NP model which contains more than two operators will also 
be allowed, as long as the measurements in both B^ (pKg and B (pK* decays 
are explained. However, models with supersymmetry or extra Z' bosons contain only 
operators of V/A type, and therefore cannot explain both B^ (pKg and B — > (pK* 
data. In contrast, the two-Higgs-doublet model has only S/P operators, and is thus 
favored. 

Acknowledgements: We thank J. Matias and M. Nagashima for helpful communi- 
cations. This work is financially supported by NSERC of Canada. 
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A Calculation of factorized matrix elements 



For NP, we work within the framework of naive factorization. Since LO contributions 
of NP are expected to be subleading compared with the SM, NLO contributions of NP 
can be safely neglected. Then, for a general effective four-quark operator O ~ X ® 1^ 
(X, Y = S,P, V, A, T or T75), the matrix element is assumed to factorize as 



{^K\0\B)^{K\X\B) {<P\Y\0) 



(19) 



where K stands for Kg or K*. {K\X\B) is calculable using known form factors; 
(01 y |0) is calculable using the 0-meson decay constants. 
For B ^ K form factors, we use definitions from Refs. 



{K{p)\V B{p, 



m 



K 



+ 



K 



(K(p)|T(l±75) B{p 



B)/ q 



'B)i q 



{K*{p,e)\V±A B{p,) 



{K{p)\T B{p, 
i{{PB+p)t,s-q^{ml -m^)} 

Ms) 



fAs) 



TJig + rriK 



tKpb +p)A^* -Pl 



Tiq^^ -Pi 
(ir*(p,e)|T(l±75) B{p,)) j' 



(A3(s)-Ao(s)) + e^,,.e*XP'^ 



2V{s) 



nig + m 



K* 



±^2(5) [el{ml-ml,) - (e* ■ Pb){Pb + p)t,} 



±Ts{s)ie*-p,)\q^- 



m'i — m?,. 



(Pb + p) 



n ( 1 



(20) 



with q = Pb — P and s = q^. Values of the form factors are tabulated in Tables 3, 4 
and 5 of Ref. [30] with minimal and maximal values and s dependance. Using these, 
we have calculated them for the case of s = (see Table [9]). 
The vector-meson decay constants are defined by [301 EI] 



(0(g,e)|y|Or 
e)|T|0)''" 



which imply 



(0(g,e)|r75|O) 



flU 



(21) 
(22) 
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Minimal value 


Central value 


Maximal value 


J + 


0.295 


0.337 


0.391 


fo 


0.286 


0.327 


0.379 


fr 


0.319 


0.375 


0.446 


Ai 


0.301 


0.345 


0.393 


A2 


0.258 


0.295 


0.333 


Ao 


0.437 


0.498 


0.750 


V 


0.423 


0.483 


0.579 


Ti 


0.355 


0.402 


0.463 


T2 


0.341 


0.387 


0.445 


T3 


0.245 


0.272 


0.307 



Table 9: Values of B ^ K, B K* form factors for s = following Ref. [30 
(calculated in the QCD light-cone sum- rules approach at the scale = rrib). 



Values of the decay constants are tabulated in Table [TUJ 



M 


[MeV] 


B [MeV] 


K [MeV] 


K* [MeV] 


f M 


215 ±5 


200 ± 25 


160 


220 ± 5 


J M 


186 ±9 






185 ± 10 



Table 10: Values of decay constants for mesons. [31j 

In order to calculate factorized matrix elements, we define the four-momemta 
Ps = (ms, 0,0,0) , pif^ji* = {E,i^ji*, 0,0, -pc) , p^ = {E^,0,0,pc) , (23) 
and polarization 4- vectors 



-{p„0,0,-E, 



-{p„0,0,E^) , 



1 

1 



(0,^1, +2,0) 



(0, +1,-2,0) , 



(24) 



in the 5-meson rest frame. From here it is straightforward to calculate factorized 
matrix elements using the above. Those which are non-zero are 



{K\V B) {<P\V\Or = 2UUm,p,, 

(J, 



for 50 ^ (PKs and 



{K*\T 



B 



111/ 



{K\T 



T\0) 



B 



liv 



T\0) 



^frfi 



"mg + rriK 



(25) 



flU 



(7^*1 T75 
15 



B 



T75 |0) 



fii/ 



A=± 



{K*\T B)^^ (0IT75IO) 



{K*\T B) (0IT75IO) 



{K*\V b) (0|1^|O)' 



{K*\A B) i(p\VW 



i,K*\A b) {(pivioy 



A=0 



A=± 



A=± 



A=0 



fi 



{-T2(m| -ml- ml,){m' 

•)!} 



+4:plml T2 + T^ml/ (ml - m 
m^rngpc 



T2tUV 



tUb + rriK* 



ttIb + rrif 



m\-m\,)Ai 



^m\p\A2 



A=± 



T B ^ (j)K* . All of this allows us to calculate the entries in Table [21 
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